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Abstract 

Three-dimensional transonic viscous flow compu- 
tations are presented for a generic high-speed acceler- 
ator model which includes wing, body, fillets, and a 
no-flow through engine nacelle. Solutions are obtained 
from an algorithm for the compressible Navier-Stokes 
equations which incorporates an upwind-biased, flux- 
vector-splitting approach along with longitudinally- 
patched grids. Results are presented for fully turbulent 
flow assumptions and include correlations with wind 
tunnel data. A good quantitative agreement for the 
forebody surface pressure distribution is achieved be- 
tween computations and the available wind-tunnel mea- 
surements at Moo = 0.9. Furthermore, it is demon- 
strated that the flow is stagnating around the boattail 
region due to separation from the aft-engine cowl lip. 

Nomenclature 

a speed of sound 

Cl lift coefficient, Lift/ tfoo^Ve/ 

C p pressure coefficient, (p — Poo)/<?oo 

Cp critical pressure coefficient 

Cp specific heat at constant pressure 

c v specific heat at constant volume 

t total body length 

E 0 total energy per unit volume 

P, (5, fi flux vectors 

h 0 total enthalpy, (E 0 + p)/p 

J Jacobian of the coordinate transformation 

k conductivity, c p pf Pr 

Moo freestream Mach number 

Pr Prandtl number 

0 state vector, J -1 ^, pu, pv, pty, E 0 ] T 

q total velocity 

qoo freestream dynamic pressure 

Rt Reynolds number based on i 
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S ref 

area of reference wing planform 

U, V, W 

body-axis Cartesian velocity components 

V* 

wall-friction velocity, \Jt w j p 

xjt 

normalized longitudinal distance aft of nose 

y + 

inner-law variable, yv* jv 

a 

angle of attack, degrees 

ft, 6 

residual parameters 

1 

ratio of specific heats, c p jc v 

P 

viscosity 

V 

kinematic viscosity, pjp 

9 

azimuthal angle 


body-fitted coordinates 

p 

density 

r w 

wall shear stress 


Introduction 

The proposed National Aero-Space Plane (NASP) 
has revived interest in hypersonic flow research in the 
past few years, particularly, in the area of computer 
algorithm developments. The design of such a vehicle, 
inevitably, would rely heavily on Computational Fluid 
Dynamics (CFD), since the current ground test facil- 
ities are limited to lower speed regimes. Also, CFD 
results generally provide much more information about 
the flow structure and its detail mechanism than one 
can obtain experimentally. 

The challenges associated with grid generation and 
surface definition along with a relatively large com- 
putational memory requirements hampers the appli- 
cations of advanced CFD-code methodologies to com- 
plex aircraft configurations. None the less, recent 
progress has been shown for several cases including 
subsonic inviscid flow computations about the com- 
plete F-14 configuration 1 as well as transonic viscous 
flow about the F-16A 2 . In addition, viscous flow com- 
putations about the ascent configuration of the space 
shuttle 3 have shown good agreement with flight data 
at subsonic, transonic, and supersonic speeds. Most 
recently, three-dimensional viscous flow computations 
about the F/A-18 forebody-LEX configuration 4 have 



demonstrated good correlation with the experimental 
wind-tunnel data as well as flight test results. Prior 
computations for generic NASP type configurations 
have demonstrated the applicability of Navier-Stokes 
methodology to transonic 6 as well as hypersonic 6 flow 
with Mach numbers exceeding 20. 

The present investigation is directed toward ap- 
plying an extended version of an implicit Navier- 
Stokes algorithm 7-10 to a generic high-speed acceler- 
ator model at — 0.9, R( « 30 X 10 6 , a *= 2°. 
At these conditions the flow is assumed to be fully 
turbulent, hence only turbulent computations are per- 
formed. The algorithm permits longitudinally-blocked 
grids which are necessary for accurately modeling the 
subject configuration. The accelerator geometry has 
been selected for this study primarily due to the avail- 
ability of current transonic wind-tunnel data 11,12 as 
well as a growing interests at NASA in utilizing the 
state-of-the-art CFD codes to analyze the flow about 
the NASP type configuration in all speed regimes. 

Governing Equations 

The governing equations as well as computational 
method for the present investigation have been pub- 
lished many times in the open literature 7-10 as they 
have evolved. The flow is presumed to be governed by 
the unsteady Reynolds-averaged Navier-Stokes equa- 
tions which are written in a body-fitted coordinate sys- 
tem. They are written in a usual conservation-law form 
as 

<5,i + (P - /.).< + (<? - <X),„ + (6 - £„),< = o 

Here the subscripts with a comma denote partial differ- 
entiation, the subscript v identifies the viscous terms, 
and the superscript A indicates scaling with respect to 
the Jacobian of the coordinate transformation. With 
the ideal gas assumption, the pressure and total en- 
thalpy can be expressed as 

p = (7-l)(£ 0 -ipg 2 ) ; 

For the present study, the thin layer approximation to 
the governing equations is invoked (i.e., P v — & v = 0) 
thus accounting for viscous flux terms only in the £ 
direction (normal to the body). Detail for these terms 
are included in Appendix 1. 

Turbulence effects are accounted for through the 
notion of an eddy viscosity and eddy conductivity: 

P - Hi + Mt - IH (1 + Ht/ni) 
k = k t + kt = k t (1 + kt/ke) = ^1 + ~p~) 


Here the subscripts i and t denote laminar and tur- 
bulent, respectively. The algebraic turbulence model 
developed by Baldwin and Lomax 13 is used to evaluate 
appropriate turbulence quantities. For separated flow 
regions, the notions of Degani and Schiff 14 are drawn 
upon to determine proper turbulence length scales. 

Computational Method 

Discretization of the governing equations results in 
a consistent approximation to the conservation laws in 
integral form 

If J Jqw + J J?hds = o 

where the time rate of change of the state vector Q 

within a cell is balanced by the net flux / across the cell 
surface. The convective and the pressure terms arising 
in the flux quantities are represented using an upwind 
biased algorithm, with either the flux-vector-splitting 
approach of van Leer 16 or the flux-difference-splitting 
approach of Roe 16 . Both approaches rely on a recon- 
struction of primitive variable data at the cell-centers 
to the cell-interfaces; the reconstruction is monotone to 
prevent oscillations near discontinuities, such as shocks, 
and is third- order accurate in the special case of one- 
dimensional flows. The shear stress and the heat trans- 
fer terms are centerally differenced so that the resulting 
algorithm is spatially second-order accurate. Solutions 
are advanced in time with a spatially-split approximate 
factorization method which, in general, solves a series 
of block 5x5 tridiagonal matrices to advance the solu- 
tion to a new time level. For the flux-difference splitting 
scheme, the implicit terms are diagonalized, leading to 
a series of scalar tridiagonal inversions. 

In the present computations, attempts were made 
initially to use flux-difference splitting with third-order 
accuracy. However, due to convergence difficulties the 
alternative flux-vector-splitting approach was invoked. 
Prior results 17 have shown the van Leer scheme to dif- 
fuse boundary layers, particularly when utilized with 
first-order accuracy and/or with sparse grid resolution. 
This limitation was overcome in the present work by 
utilizing the van Leer scheme with third-order accuracy 
and by performing the computations with a sufficiently 
dense mesh. Salient aspects of the van Leer formulation 
are discussed in Appendix 2. 

Information is exchanged between the longitudi- 
nal blocks using the patching algorithm described by 
Thomas, et. al 10 . At the interface between the two 
blocks, the dependent variables at the cell centers of 
each zone are interpolated across the coincident inter- 
face boundary assuming a biquadratic variation across 
each cell of the opposing zone. The interpolation is sec- 
ond order accurate and conserves the mass, momentum, 



and energy across the zonal interface only to within the 
truncation error of the discrete solution. 

Experimental Models 

There were two accelerator models tested in the 
Langley’s 16-foofc transonic wind-tunnel. The first 
model was a sting mounted model with an open in- 
let and exhaust faces (i.e., flow-through). This model, 
called the “forebody” model, was intended to only pro- 
vide experimental pressure data on the forebody. The 
second model, called the “propulsion” model, had a 
faired over (closed) inlet and open exhaust face to al- 
low for the simulation of propulsion effects. This model 
was mounted on a strut which housed some of the ex- 
haust jet-flow propulsion apparatus. The experiment 
was designed to obtain experimental pressure measure- 
ments on the boattail region for various exit mass-fluxes 
including zero (i.e., a jet-off case). 

Numerical Model Representation 

The surface grid definition is obtained from an an- 
alytical representation of the configuration. The com- 
plete surface representation of the configuration along 
with various longitudinal cuts are shown in Fig. 1. The 
body of the configuration consists of a 5° half-angle 
right-circular cone and a 9° boattail frustum that are 
connected with a cylinder. A nacelle is wrapped around 
the cylindrical part of the body. The configuration in- 
corporates wing-fillets on both the starboard as well as 
the port side of the body. The front engine wing-fillet 
blends into the nacelle cowl at a longitudinal station be- 
tween stations C and D (Fig. 1). The aft-engine wing- 
fillet merges smoothly into the wing upper and lower 
surfaces, at a longitudinal station between stations E 
and F. The configuration included a 70° swept sharp 
delta wing with a span of 0.11 X t. The wing had a 4% 
thick biconvex airfoil and was mounted at 3° incidence 
on the combinations of front and aft-engine fillets as 
well as the nacelle. 

The inlet face is faired over similar to the tested 
wind-tunnel propulsion model. Since the effects of 
propulsion are not included in the present computa- 
tions, the exhaust face was also faired over with a 
45° ramp. This fairing eliminates the problem of hav- 
ing a surface discontinuity at the exit plane, where 
there needs to be a longitudinal patching station of 
two blocks from a geometrical consideration (i.e., start 
of the aft-fillet). The surface geometry is defined with 
a total of 5,400 grid points at 88 longitudinal stations 
along the body. 

The field grid is divided into eight longitudinal 
blocks, each representing a local geometrical complex- 
ity of the configuration. The forebody block is gener- 
ated with a C-0 grid topology thus resolving the nose 


radius and extending to the upstream inflow boundary. 
Figure 2 shows the forebody grids on the surface, plane 
of symmetry, as well as the base section. In this figure, 
the radial extent of the grid has been truncated for 
clarity. The three-dimensional grid for the remaining 
blocks are constructed from two-dimensional O-type 
cross flow grids which are longitudinally stacked, con- 
stituting an H-0 topology. Figure 3 illustrates a farfield 
as well as a nearfield view of complete grid about the ac- 
celerator configuration in the plane of symmetry. The 
grid lines in the radial direction are plotted for every 
fourth point for clarity. The grid dimensions as well 
as the mesh topology for each corresponding block are 
listed in Table 1. 


Block 

X 


r 

Points 

Top 

1 

29 

31 

65 

58435 

C-0 

2 

12 

59 

65 

46020 

H-0 

3 

12 

91 

65 

70980 

H-O 

4 

9 

79 

65 

46215 

H-0 

5 

13 

91 

65 

76895 

H-O 

6 

9 

75 

65 

43875 

H-O 

7 

4 

31 

65 

8060 

H-O 

8 

11 

31 

65 

22165 

H-0 


Table 1. Grid dimensions and topology. 


The upstream, downstream, and the radial extent 
of the flowfield grid is about one configuration body 
length which corresponds to 13.5 times the maximum 
body diameter. Block 8 is the downstream extention 
of the grid and has 11 longitudinal stations. All blocks 
have 65 points in the radial direction and the circum- 
ferential number of grid points vary from block to block 
depending on the complexity of the local cross-sectional 
geometry. Longitudinally, the grids are clustered near 
each geometrical break such as a station where the wing 
starts or terminates. There axe approximately 373,000 
grid points used to represent the entire flowfield do- 
main. 

The crossflow grids are generated using established 
transfinite interpolation techniques 18 ’ 19 with a method 
applicable to slender shapes 4,6 . The flowfield grid is 
generated with sufficient normal clustering near the 
surface to adequately resolve the laminar sublayer of 
the turbulent boundary layer flow at the subject wind- 
tunnel freestream conditions (Moo = 0.9, Re fa 30xl0 6 , 
a — 2°). This grid produced an average normal cell- 
center size next to the wall of approximately 1 0~ & £ 
which corresponded to y + » 3 for the turbulent com- 
putations; a laminar sublayer generally extends out to 
y + fa 8.5. Fig. 4 illustrates the nearfield view of a typ- 
ical configuration cross section as well as the grid res- 
olution in the vicinity of the fillet-body and fillet- wing 




juncture. This figure illustrates the challenge that is 
associated with generating a single O-type grid around 
such complex cross section with various break points in 
the surface geometry. 

Results and Discussion 

Results are presented first for general flow features. 
This is followed with a discussion of some convergence 
attributes of the computed results. Finally, a compar- 
ison between the computed surface pressures arid the 
available experimental data are presented on the fore- 
body as well as boattail region. 

Solution Attributes 

Mach contours on the surface and in the config- 
uration plane of symmetry at the subject wind-tunnel 
flow conditions are shown in Fig. 5. The sonic line 
is highlighted with a white contour line in the plane 
of symmetry to highlight the supersonic flow region. 
The magnitudes associated with contour quantities are 
displayed with a color bar. This figure indicates that 
after a basically subsonic forebody, the flow accelerates 
supersonically over the faired-over inlet, follows with 
a mild compression on the engine cowl and it subse- 
quently shocks down at the aft-engine cowl-lip. The 
presence of the shock along with the slanted exhaust- 
face produces a large adverse pressure gradient which 
causes the flow to separate at the exhaust cowl-lip and 
subsequently envelop the entire boattail region. Note 
that the Mach contours smoothly cross over the many 
block interfaces. Cross-flow Mach contours are shown 
in Fig. 6 for various planes along the length of the con- 
figuration. This figure illustrates the radial extent of 
the supersonic zone in two cross-flow planes: one on the 
nacelle and the other slightly aft of it cutting through 
the shock. In addition, the complex flow structure in 
the separated boattail region is well depicted in the 
last two cross-flow planes. It should be noted that the 
propulsion effects (i.e., either a simple flow-through or 
the actual engine simulation) are expected to have a 
significant influence in the supersonic region as well as 
the separated flow in the boattail region. 

The total pressure contours in various cross-flow 
planes along the body as well as in the plane of sym- 
metry are shown in Fig. 7(a). The flow appears to be 
primarily attached on the forebody. Also, the losses 
that are associated with viscous flow are well depicted 
in the separated boattail region. Figure 7(b) shows a 
closeup view of the total pressure contours to highlight 
the smaller separated flow structures that exist on the 
engine cowl as well as the wing upper surface. The sep- 
arated flow on the cowl is associated with the fillet-body 
juncture flow that is spilling over the cowl on its lower 
surface and, to a lesser degree, on its upper surface. 


The small flow structure on the wing upper surface 
is evidence of a leading-edge bubble type separation. 

The flow intensity at the center of this separated flow 
is reduced significantly as it passes through the shock ^ 
which is sitting at the aft-engine cowl-lip and extended 
circumferentially over onto the wing surface. 

The computed streamlines on the configuration 
upper surface as well as in the plane of symmetry are 
presented in Fig. 8. This figure clearly shows the abrupt 
change in the streamline pattern on the wing upper 
surface near the aft-engine cowl-lip that signifies the 
presence of the shock. The turning of the surface flow 
streamlines is primarily due to the diminished axial ve- 
locity component of the flow. Consequently, in the 
wing tip region, the stream fines appear to go around 
the wing leading edge to the lower surface. Fig. 8 also 
shows the very complex surface flow pattern that is as- 
sociated with the low speed flow in the boattail region; 
this pattern appears to indicate a whorl. 


The flow computations were performed on the 
NAS Cray 2, located at NASA-Ames. On this ma- 
chine, the algorithm requires approximately 50 fx sec- 
onds per grid point per cycle. The present numerical 
results for the a — 2° were obtained in nominally 4815 
cycles which required about 25 hours of computer time. 
This number of cycles was sufficient to reduce the resid- 
uals by three orders of magnitude and limit the oscil- 
lations in Cl, to « ±.005 (Fig. 9). The computations 
were performed without the use of mesh sequencing or 
multigrid iteration. 

The salient features of the flow were examined dur- 
ing the course of the convergence process. It was found 
that the position of the shock and the shape of the su- 
personic pocket was apparently established within ap- 
proximately the first 2000 cycles. The loss in Cl, shown 
in Fig. 9, for 1500 < iterations < 2000 is primarily as- 
sociated with the formation of separated flow in the 
boattail region. Subsequent to establishing the shock 
structure (iteration « 2000), changes in this separated 
flow region result in a further loss of lift for the next 
ps 1000 cycles after which changes in lift become small. 
The separated flow region appears to exhibit certain 
unsteady flow characteristics which axe associated with 
the downstream convection of flow quantities from the 
exhaust face. There are reasons to believe that this sep- 
arated region is the primary cause of the oscillations in 
the lift as well as the residual magnitude for the last 
1500 cycles (see Fig. 9). 

To examine this conjecture, the individual block 
rms residuals ($„) tire appropriately weighted by the 


Convergence Characteristics 


4 



fraction of grid points in each block, 




npn 

. n Ptotai , 




where 


ifin) = 


*=np n 


1 v-A 

) residual ? 

n Pn fr[ 


In the above equations, residua^ refers to the resid- 
ual in cell i and np„ refers to the number of points 
for block n. The total rma residual for the complete 
configuration can now be obtained from 


residualtotoi — 


\ 


«— i 


where N is the total number of blocks. The block resid- 
ual distribution is presented in Fig. 10 for four different 
levels of iterations. This figure shows that the incre- 
mental local residual for blocks 5, 6, and 7 are rela- 
tively more dominant than the corresponding residuals 
for the upstream blocks 1 through 4. The separation in 
the boattail region is mainly responsible for the large 
local residuals as well as their oscillatory behavior with 
number of iterations. This separation also appears to 
be the principal cause for the excessive total number 
of iterations performed to achieve the present results. 
Note that the local block-residual for the first 4-blocks 
are not changing much by advancing the solution from 
2455 to 4815 cycles; the flow on the forebody appears 
to have setup within the first 2455 cycles. 

The surface pressure sensitivity on the develop- 
ping solution was also examined at three iteration lev- 
els. Figure 11 illustrates longitudinal pressure distribu- 
tions on the accelerator forebody at five conical rays; 
0=0° corresponds to the windward and 6 = 180° to 
the leeward side of the configuration plane of symme- 
try. This figure also shows a schematic side view of the 
computational accelerator body to highlight the longi- 
tudinal extent of the forebody. As discussed previously 
in the solution attributes, the forebody appears to ex- 
hibit a fairly benign flow condition with the expected 
flow compression in the vicinity of the inlet face. It is 
also evident that the computed surface pressures are 
less sensitive to the last 1310 cycles beyond 3505. 

The boattail surface pressure sensitivity in the 
plane of symmetry for the same number of cycles are 
shown in Fig. 12(a). This is basically the continua- 
tion of the forebody pressures shown in the previous 
figure onto the engine-cowl and the boattail region in 
the plane of symmetry. The flow accelerates over the 
faired-over inlet supersonically (C* (a —0.2) and then it 


goes through a mild compression over the engine cowl 
and subsequently shocks down at the aft- engine cowl- 
lip. Note that the supersonic region, as well as the 
shock strength and its position are well established with 
r* 2800 cycles. However, the pressures on the boattail 
region appear to be somewhat sensitive to those num- 
ber of cycles examined. 

The circumferential pressure distribution at two 
stations, one ahead of the shock (x/l = 0.768) and the 
other aft of the shock (x/l = 0.828), are presented in 
Fig. 12(b), along with their corresponding geometrical 
cross-section. The left part of the figure shows the vari- 
ation of pressures as a function of azimuthal angle $ on 
the body. The right part of the figure shows the pres- 
sures on the wing ( x/l = 0,768) and wing as well as the 
fillets (x/l = 0.828) as a function of the exposed wing 
or wing-fillet semispan, respectively. Notice these plots 
have scales that are three times as fine as the figures on 
the left. The pressure distribution at the station ahead 
of the shock appears to be insensitive to the number of 
cycles, particularly on the wing. Furthermore, the ef- 
fects of the separated flow emitted from the fillet-body 
juncture over the engine cowl surface (x/l = 0.768) are 
well evident at 6 w 60° (i.e. lower surface) and 9 « 100° 
(i.e. upper surface). As discussed earlier, the pressures 
at the station aft of the shock experiences some sen- 
sitivity to the number of cycles. The pressures on the 
wing clearly show the loss of the wing leading-edge suc- 
tion after passing through the shock. 

In summary, the results from the surface pressure 
sensitivity study indicate that the flow on the forebody 
appears to had been well established within the first 
3500 cycles. These results also reveal that the super- 
sonic pocket over the engine cowl as well as the shock 
position and its strength were fully developed at about 
2800 iterations. Furthermore, the flow in the separated 
boattail region were most sensitive to the number of 
cycles. Note that these results complement the find- 
ings discussed earlier from the residual-history point of 
view. 

Data Comparison 

The computed static pressure coefficients on the 
forebody are compared with the experimentally mea- 
sured data 11,12 in Fig. 13, The figure shows the com- 
parison at five different conical rays from windward 
(9 = 0°) to leeward side (9 = 180°) in 45° increments. 
A schematic side view of the accelerator body for both 
the wind-tunnel model as well as the computational 
model are also shown to highlight the major geometri- 
cal differences that exist between them. As described 
previously, the experimental fore body pressures were 
obtained on the same model but only for flow through 
inlet conditions whereas the numerical model had the 
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faired over inlet and exhaust face. The longitudinal ex- 
tent of the forebody covers approximately 70% of the 
configuration total body length. This figure reveals a 
good correlation between experimental data and the 
computed results up to x(t » 0.4. The discrepancies 
between theory and data in x/i range of 0.4 - 0.7 are 
mainly due to the differences in inlet flow conditions. 
The faired over inlet produces additional compression 
which results in a more positive pressures as indicated 
in the computed results of Fig. 13. 

The computed pressure coefficient on the boattail 
are plotted against the measured data in Fig. 14(a), 
along the configuration plane of symmetry. This figure 
also shows a schematic side view of the experimental 
as well as computational model to illustrate the major 
differences that exist in the exhaust face with regard to 
the geometry and the flow condition. The data corre- 
spond to the experimental condition of zero-mass flux 
across the exit. This figure reveals a fairly good corre- 
lation between theory and experiment. The discrepan- 
cies in the neighborhood of the exit could well be due 
to geometrical differences at the exit face and/or the 
back pressure effects from the open exhaust face in the 
experiment. 

The circumferential theoretical and the experimen- 
tal data for two longitudinal stations on the boattail 
region are shown in Fig. 14(b). The first station is lo- 
cated slightly aft of the shock [xfl = 0,828) and the 
other further downstream [xjl — 0.961). Notice that 
the scales have been magnified to bring out the differ- 
ences. At xjl — 0.828 there are some disagreements 
between theory and experiment although the trend ap- 
pears to have been better predicted. At xfl — 0.961 
the agreement between theory and experiment is quite 
good. In general, the flow structure in the boattail 
region is very complex, and the correlations that have 
been achieved there are perhaps better than might have 
been expected. 


Concluding Remarks 

Three-dimensional transonic viscous flow compu- 
tations for a complex high-speed accelerator are pre- 
sented with fully turbulent flow assumptions. Good 
correlation between experimental static pressures and 
computational results on the forebody are disclosed. 
The computed results show a large supersonic pocket 
over the nacelle surface with a shock at the aft-engine 
cowl-lip. The presence of the shock as well as the 
slanted exhaust face cause the flow to separate into 
a very complex flow structure envelopping the entire 
boattail region. None the less, some quantitative com- 
parisons have been made between theory and experi- 
ment in the boattail region. 


Appendix 1 - Governing Equations 

The inviscid flux terms from the governing equa- 
tion are defined as follows: 


' p ' 


pU \ 

pu 

. 1 

pUu + 

pv 

: ^=4 

pUv + ivP 

pw 

J 

pUw + ^,p 

> E a , 


{ ( E 0 +p)U J 


pV 


( pw 

pVu + r}, x p 

. i 

pWu + 

pVv+ r}, y p 

; * = - 

pWv + $ <t) p 

pVw-\-r}' Z p 

J 

pw w + $,zp 

( E 0 +p)V J 


{ {E 0 +p)W 


The contravariant velocity components in these 
equations are given by: 


u ) ( S,* u + tv v + 1 

V > = < »?,*« + V, v v + V,*u> f 

W ) [ f ti x u + $ iV v + f t ,u> j 


The thin layer viscous flux terms from the govern- 
ing equation are defined as follows: 


A,- 


IhT 


0 

i + s,x4>{ 2 

< + fojfg ► 

2 

. (( V).f + Prf-1- + Wfa , 




where 


{*«i> = {(s.) a + fo.) 2 + fe*n 


{ 4(1 } = 5 { “.(O + } 

The transformation metrics are: 


V,x $,x 
£,y V,y 

L J 


- iU*,* ~ 

z ,*t x ,s — x ,*} z ,s x ,i z ,< ~ z ,t x ,$ z ,Z x ,f) ~~ x ,t z >i) 
L x ,nV,f ~ y,n x ,s y>z x ,( ~ x ,zy,s x ,zy,n ~ y,t x ,n 

and the Jacobian of the transformation is given by: 

J ~ l = *Av,n*, t - ~ - *,**,*) 

+ z A x >ny,< - y,n x ,<) 
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Appendix 2 - Upwind Algorithm 
Generalized convective and pressure fluxes are split 
into forward and backward contributions and differ- 
enced according to the van Leer flux-vector splitting 
scheme. Using the £— direction as an example 

{sf'/Hh = [f(P + +s^-]i 

= [/+(«-) + #-(<? + ) 

-i^ + (<n i b 

where — 1 and ^ + (Q - ) <+ denotes a forward 
flux evaluated with an upwind biased interpolated state 
variables. 

The flux is split according to the contravariant 
Mach number, Afj = e/a, where 

o s 17/ | yrad(£) | 


with 

minmod(x f y) s 3 tyn(x) x 

max [0,rm'n(x sign(y),y siyn(x))] 

Here B is a parameter which controls the maximum 
allowable range of successive gradients before limiting 
is applied. It is defined in terms of k 

B ~ (3 — /c)/(l — /c) 

and k = 1 /s for third-order accuracy in the special case 
of one-dimensional flow. 
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For supersonic flow | \ > 1 and 

P + = P } p~ = 0; M i > +1 

P~ = P, P + = 0 ; M t <- 1 


For subsonic flow | j< 1 and 


/ ± 

p± _ J mass 


|yrad(e) 


where 


|,*(~0±2a)/7 + u 
| iV (~0±2a)/7 + u 
dz 2a)/7 + w 

f± 

J energy 


fL.. = ±P°(M t ± 1)74 


ftnerm = [-(l-l)o 2 ±2('y-l)na+2a 2 )/(i 2 -l)+s72 
and 


\lx, iy , f,,] T = \ix, ivy U T / I grad[t) I 

State variables at the cell interfaces are formed 
from interpolation of the primitive variables vector 
q = (p, u, v, w,p] T by the one-parameter family 

(?~)»+i/2 = 7* + ‘AH 1 “ *)V7» + (1 + /c ) A 9t] 

h + )<+ Vi " qi + l ~ V^K 1 + k)V9»+i + (! “ 

where the backward and forward differences are ex- 
pressed respectively with a typical minmod limiter as 

V? — rmnmod[\/q } HA g] 

A q — minmod\Aq , 
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Figure 2.- Accelerator forebody C-O grid. 





(a) Farfield. 

Figure 3.- Blocking strategy for the accelerator configuration. 
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(b) Nearfield. 
Figure 3.- Concluded. 





Figure 4.- A typical accelerator crossflow grid. 
10 






Figure 6.- Mack contours in cross-flow planes. M co =-~ 0.9, Rt 30 X 10 6 , a ~ 2°. 
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Figure 10.- Block residual history. Moo = 0.9, R& R* 30 x 10 6 , a ~ 2°. 
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Figure 12.- Boat tail surface pressure sensitivity on solution development. 
Moo = 0.9, R t xi 30 x 10°, a = 2°. 
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Figure 13.- Forebody surface pressure correlations with experiment. 
Moo ~ 0.9, R t « 30 x 10 6 , a — 2°. 
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